Charge-imbalance effects in intrinsic Josephson systems 



S. Rother, Y. Koval, P. Miiller 
Physikalisches Institut III, Universitat Erlangen-Niirnberg, Erlangen D-91058, Germany 



(N 
O 
O 



R. Kleiner 

Physikalisches Institut, Universitat Tubingen, D-72076 Tubingen, Germany 

D. A. Ryndyk, J. Keller 
Institut fur Theoretische Physik, Universitat Regensburg, D-93040 Regensburg, Germany 



(N 

O 
O 

Or 



C. Helm 

Institut fur Theoretische Physik, ETH Zurich, CH-8093 Zurich, Switzerland 
(Dated: February 1, 2008) 

We report on two types of experiments with intrinsic Josephson systems made from layered high- 
T c superconductors which show clear evidence of nonequilibirum effects: 1. In 2-point measurements 
of IV-curves in the presence of high-frequency radiation a shift of the voltage of Shapiro steps 
from the canonical value V s = hf/(2e) has been observed. 2. In the IV-curves of double-mesa 
structures an influence of the current through one mesa on the voltage measured on the other 
mesa is detected. Both effects can be explained by charge-imbalance on the superconducting layers 
produced by the quasi-particle current, and can be described successfully by a recently developed 
theory of nonequilibrium effects in intrinsic Josephson systems. 

PACS numbers: 74.80Dm, 74.40.+k, 74.50.+r, 74.72.Fq, 74.72.Hs 
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I. INTRODUCTION 

In the strongly anisotropic cuprate superconductors 
Bi 2 Sr 2 CaCu20 8+( 5 (BSCCO) and Tl 2 Ba 2 Ca 2 Cu 3 Oi 0+ 5 
(TBCCO) the Cu0 2 layers together with the interme- 
diate material form a stack of Josephson junctions. In 
the presence of a bias current perpendicular to the lay- 
ers each junction of the stack is cither in the resistive or 
in the superconducting state leading taj-thc well-known 
multibranch structure of the IV-curves Buu 

In the case of weakly coupled layers we expect that the 
bias current generates charge accumulation on the layers 
between a resistive and superconducting junction. As the 
current through a resistive junction is carried mostly by 
quasi-particles, while the current through a barrrier in 
the superconducting state is carried by Cooper-pairs, we 
expect two types of non-equilibrium effects: 1. charge 
fluctuations of the superconducting condensate, which 
can be expressed by a shift of the chemical potential 
of the condensate and 2. charge-imbalance between 
electron- and hole- like quasiparticles. 

Non-equilibrium effects in layered superconduc- 
tors have been discussed in a number of pa- 
pers in variouSpfiOJitfYKl.s ..with different methods and 
approppiations.QHD'OQBEjO In a recent theoretical 
papertL2l we have investigated in particular the conse- 
quences of non-equilibrium effects on experiments with 
stationary currents. We found that in this case only 
the charge-imbalance is important and leads to a change 
in the voltage, while the shift of the chemical potential 
has no influence on the measured voltage. The latter 
determines e.g. the dispersion of longitudinal Joseph- 
son plasma wavesu and can be observed in some optical 
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In this paper we report on new experiments which show 
clear evidence of non-equilibrium effects for stationary 
currents and which can be explained by charge-imbalance 
on the superconducting layers. In the first type of ex- 
periments Shapiro steps produced by high-frequency ir- 
radiation are measured in mesa structures of BSCCO 
with gold contacts. Here a shift SV of the step-voltage 
AVs = hf I (2e) — SV from its canonical value hf / (2e) is 
observed, which can be traced back to a change of the 
contact resistance due to charge-imbalance on the first 
superconducting layer. In another type of experiments 
current- voltage curves are measured for two mesas struc- 
tured close to each other on the same base crystal. Here 
an influence of the current through one mesa on the volt- 
age drop on the other mesa has been measured which can 
be explained by charge-imbalance on the first common 
superconducting layer of the base crystal. Both experi- 
ments allow to measure the charge-imbalance relaxation 
rate. .— . 

We start with a brief summary of the theoryta which 
will be used in the following. Then the sample prepa- 
ration and the different experiments are described and 
discussed. From the experiments the charge-imbalance 
relaxation time will be determined. 



II. OUTLINE OF THE THEORY 

Let us consider a stack of superconducting layers n = 
1,2... with a normal electrode n = on top. The ba- 
sic quantity which determines the Josephson effect of a 
junction between layer n and n + 1 is the gauge invariant 
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phase difference 

2e f n+1 

7n,n+l(i) = Xn(t) - Xn+lW ~ J dzA z (z,t), (f) 

where Xn(t) is the phase of the order parameter on layer n 
and A z (z, t) is the vector potential in the barrier. For the 
time derivative of J n ,n+t one obtains the general Joseph- 
son relation: 

2e / \ 

7n,n+l = Y (K,n+1 + $„+i - $„J . (2) 

Here 

rn+l 

V n ,n+1= dzE z (z,t), (3) 

•/ n 



= " ^-Xn(*), (4) 

2e 

are the voltage and the gauge invariant scalar potential, 
4> n {t) is the electrical scalar potential. The quantity e<&„ 
can be considered as shift of the chemical potential of the 
superconducting condensate (in this paper the charge of 
the electron is written as — e). 

The total charge fluctuation 8p n on layer n consists 
of charge fluctuations of the condensate and of charge 
fluctuations of quasi-particles. It is convenient to express 
the latter also by a kind of potential writing 

Sp n = -2e 2 iV(0)($„-* n ). (5) 

With help of the Maxwell equation (d is the distance be- 
tween the layers, e the dielectric constant of the junction) 

fiPn = -^J-{V n . n +l ~ V n -l, n ) (6) 

a 

the generalized Josephson relation now reads: 

h , 

2e 

+ (7) 

with a = eeo/ ' (2e 2 N(0)d). It shows that the Joseph- 
son oscillation frequency is determined not only by the 
voltage in the same junction but also by the voltages in 
neighboring junctions. Furthermore it is influenced by 
the quasi-particle potential ^ on the layers. If we ne- 
glect the latter we obtain for j n ,n+i the same result as 
in Ref. | 

These equations for the voltage between the layers have 
to be supplemented by an equation for the current den- 
sity: in the stationary state (no displacement current) it 
can be written asE3 

jn,n+i — j c sin7„ in+ i 
+ ^ b f ± ((1 + 2a)V n<n+ i - a{V n -i, n + V n+ i, n+2 )) . (8) 



Here the quasi-particle current between layers n and n+1 
is driven not only by the voltage V n>n ±i between the lay- 
ers but also by additonal terms which result from the 
charge fluctuation on the two layers. In the stationary 
state the dc-dcnsity j n ,n+i is the same for all barriers 
and is equal to the bias current density j. If the Joseph- 
son junction is in the resistive state we may neglect the 
dc-component of the supercurrent density in (^) for junc- 
tions with a large McCumber parameter j3 c 3> 1. Then 
we obtain: 

jd 

= (1 + 2a)V n . n+ i - a(V n -l,n + K+l.n+2)- (9) 

0"n,n+l 

For the junction between the normal electrode (n = 
0) and the first superconducting layer (n — 1) we may 
neglect charge fluctuations on the former and obtain for 
the current: 

il = (1 + a)V ,i - aVi, 3 . (10) 
fo,i 

Finally we need an equation of motion for the quasi- 
particle charge. Here we consider a relaxation process be- 
tween quasij-particle charge and condensate charge within 
the layer .L3E£(EZl In the stationary case the quasi-particle 
charge is proportional to the charge-imbalance relaxation 
time T q and the difference between supercurrents flowing 
in and out the layer, oij_cquivalcntly, by the difference in 
quasi-particle currents £3 

*n = C/ c sin7„ !Tl+ i - 7 c sin7 n _ Ml )T 9 /(2e 2 iV(0)) 
= {3f-x,n-iZnM{^ 2 N{Q)) (11) 

In the limit of small non-equilibrium effects and for 
T <C T c we may use the approximation j^ n+ i — 
Cn,n+iKi,n+i = j for a resistive junction and j qp ~ for 
a junction in the superconducting state. For example, if 
a current j is flowing from a junction in the resistive state 
into a junction in the superconducting state, the charge- 
imbalance potential generated on the superconducting 
layer between the two junctions is ^ n = jT q /(2e 2 N(0)). 



III. SAMPLE PREPARATION AND 
EXPERIMENTAL SET-UP 

For measuring the intrinsic Josephson effect in BSCCO 
we used a mesa geometry. Our base materials were 
BSCCO single crystals which were grown by standard 
melting techniques. After glueing these crystals on 
saphire substrates, a 100 nm thin gold layer was ther- 
mally evaporated. The mesa was patterned by electron 
beam lithography and etching by a neutral Ar atomic 
beam. As etching rates for gold and BSCCO are known 
quite exactly, the number of Josephson junctions inside 
the mesa could be adjusted to be between 5 and 10. To 
avoid shortcuts between Au leads and the superconduct- 
ing base crystal an insulating SiO layer was evaporated. 
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SiO 
Au 




Samples 


area {p,m 2 ) 


#jj 


experiments 


SR102.3 10 x 10 (B) 


6 


DM 




4 x 10 (Ml) 


10 


DM 




5 x 10 (M2) 


10 


DM 


SH104 


7 x 6 (B) 


11 


DM 




2x6 (Ml) 


3 


DM 




3x4 (M2) 


4 


DM 


#32 


8x8 


8 


FIR 


#39 


10 x 10 


20 


FIR 


#20 


10 x 10 


11 


FIR 



TABLE I: Table of samples described in this article. The area 
of the mesas and the number of junctions (#jj) containing it 
are listed together with the kind of experiments they were 
used for. 



FIG. 1: Sketch of the sample geometry used for injection 
experiments. 

This material was removed from the top of the mesa by 
liftoff. For contact leads a second Au layer was evapo- 
rated. The shape of these leads was defined by standard 
photolithography and etched by Ar atoms. For samples 
used in FIR experiments these-Au leads were formed in 
the shape of bow-tie antennas Jl3 

A second type of experiments is concerned with in- 
duced charge imbalance. Here additional preparation 
steps were necessary. These samples consist of two small 
mesas on top of one bigger mesa (Fig. [l]) . To realize this 
geometry the top layers of a larger mesa (structured as 
described above) had to be cut in two smaller ones. The 
whole sample was protected by electron-beam resist ex- 
cept a thin line on top of the big mesa. By etching this 
sample the larger mesa could be divided in two smaller 
mesas. To keep a part of the base below the two smaller 
ones the etching process had to be controlled precisely. 
The resistance of the gold lead was recorded during etch- 
ing showing the time when gold was removed and separa- 
tion of the big mesa began. With this technique we were 
able to structure samples with 5 x 10/im 2 and 4 x 10/im 2 
mesas each including about 10 intrinsic Josephson junc- 
tions on top of one larger one (10 x 10/im 2 , 6 junctions). 
The gap between the two top mesas was 1/j.m. All sam- 
ples discussed in this article are listed in Tab. ||. The 
kind of experiments they are used for are marked by DM 
(double mesa) and FIR (FIR absorption). 

IV-characteristics of our samples were recorded by ap- 
plying dc currents and recording voltages across mesas by 
digital voltmeters. If not mentioned otherwise all mea- 
surements were performed at a temperature of 4.2 K. 

As radiation source for experiments presented in Sec- 
tion 4.1 we used a far-infrared laser which was optically 
pumped by a CO2 laser. To minimize power losses we 
used a polyethylene lens producing a parallel beam. In- 
side the optical cryostat the samples were fixed in the 
center of a silicon hyperJaemispherical lens to focus the 
radiation onto the mesa.ES For some samples we used a 



very sensitive setup to detect Shapiro steps and to de- 
termine voltage of these resonances very precisely. For 
this purpose the laser beam was modulated by an opti- 
cal chopper with a fixed frequency. With the chopper 
frequency as external reference, the voltage across the 
junction was connected to the input of a lock-in amplifier 
(LI). As the LI analyzes voltage changes occurring with 
the chopper frequency the output signal Vli exhibits a 
point-symmetric structure. Thus the LI output signal 
represents the voltage difference between IV- character- 
istics with and without laser radiation. The voltage of 
the Shapiro step can easily be identified as the symmetry 
point of this signal. 



IV. EXPERIMENTAL RESULTS 
A. Shapiro steps 

Our experiments with high frequency electromagnetic 
radiation have been carried out with three samples #32 
(8 x 8fim 2 ), #39 (10 x 10/im 2 ) and #20 (10 x lO^m 2 ). 
We first discuss results for sample #32 which consists of 
eight intrinsic Josephson junctions as can be seen in the 
IV-characteristics shown in Fig. ||. 

The increasing values of critical currents may be ex- 
plained by the inhomogencous etching process during 
fabrication producing layers with increasing areas from 
top to bottom of the mesa. Then it is natural to as- 
sume that the values of critical currents increase with 
the position of the junctions inside the mesa. In partic- 
ular, the first resistive branch of the IV-characteristics 
should be assigned to the uppermost Josephson junc- 
tion of the mesa. The critical current on branch num- 
ber is strongly suppressed and its IV-curve is linear, 
while the other branches show the typical non-linear IV- 
dependence characteristic for a tunneling junction be- 
tween two superconducting layers with d-wave order pa- 
rameter. The special behaviour of branch might be 
explained by the assumption that the first superconduct- 
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FIG. 2: IV characteristic of sample #32. 



ing layer is in proximity contact with the normal gold 
electrode. 

This sample was irradiated with four external frequen- 
cies between 584 GHz and 762 GHz. On the IV charac- 
teristic we could detect first order Shapiro steps on the 
first resistive branch at an absolute voltage V$ and cur- 
rent Is- To compare the step- voltage with hf/(2e) as 
predicted by the second Josephson relation the contact 
voltage Vo{I) between Au leads and mesa in the super- 
conducting state without radiation had to be measured. 
The value AVs = Vs{Is) — Vq(Is) can then be compared 
with hf I (2e). As shown in Fig. AVs evaluated for this 
sample is strictly lower than hf/(2e) for all four frequen- 
cies. The relative shift is approximately —3%. 

This downshift can be explained if we assume that the 
Josephson junction in the resistive state which is locked 
to the external radiation with frequency / is close to 
the normal electrode, i.e. between the first and second 
superconducting layer. For this junction we have 



2e 



(1 + 2a)V h2 - a{V Q>1 + V 2j3 ) + * 2 



(121 



while for the other junctions, which are in the supercon- 
ducting state, we use 

= (1 + 2a)V„, n+ l ~ a(V n -l,n + V„ + l,„ +2 .) + *n+l - *n 

(13) 

for n > 2. Adding up these equations together with the 
current relation for the contact with the normal electrode 
( ^0| ) we obtain for the total voltage: 



V s = 



00,1 



2e' 



(14) 



Here we have assumed that ^„ = for n > 3. The contri- 
bution of ^2 drops out. Finally, the charge-imbalance po- 
tential vanishes on the first layer since the on- and off- 
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FIG. 3: Voltages of Shapiro steps measured on IV charac- 
teristic of sample #32 compared to hf/(2e) (for the higher 
frequencies the current axis is shifted). 



We have to compare this result with the voltage mea- 
sured in the absence of high-frequency irradiation, when 
all junctions are in the superconducting state. Then (13) 
holds for n > 1. Adding up now equations ( |l3| ) and (To|) 
we obtain for the total voltage (contact voltage in the 
superconducting state): 



Vo 



0o,i 



(15) 



The quasi-particle potential on the first superconducting 
layer is now given by 4"i = jT q /(2e 2 N(0)) while "t> n = 
for n > 2. Subtracting the measured contact voltage 
from the voltage of the Shapiro step we obtain for the 
step- voltage: 



AVs 



TlD 

~2e 



5V 



(16) 



flowing quasiparticle currents are equal, 



with SV = jT q /(2e 2 N(0)). The shift is proportional to 
the life-time of charge-imbalance. 

Sample #39 consists of 20 intrinsic Josephson junc- 
tions. For this sample the distribution of critical currents 
is more homogeneous (Fig. ^) than for #32 making it im- 
possible to determine the position of the junction gener- 
ating the first resistive branch inside the mesa. By using 
the sensitive measurement technique with a pulsed laser 
beam we were able to detect Shapiro steps on the first 
resistive branch at three FIR frequencies between 1.40 
THz and 1.63 THz. The voltage differences AVs of these 
resonances were calculated as described above. In con- 
trast to #32 for this sample no deviations from hf/(2e) 
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FIG. 4: IV characteristic of sample #39. 

were measured. In view of the theory presented above 
this means, that the resistive junction which is locked to 
external radiation is not close to the normal electrode or 
the charge-imbalance time is rather short for this sample. 

Finally we want to focus on sample #20 consisting of 
11 intrinsic Josephson junctions. The IV-curves of this 
sample show a very homogeneous distribution of critical 
currents making it impossible to assign any branch to a 
junction at a certain position inside the mesa. However, 
for different current cycles (increase I\ to its maximum 
value than decrease it zero) two different first branches 
marked la and lb in Fig. @ are measured. This means 
that two different junctions in the stack become resistive 
first. 

We were able to detect Shapiro steps at five frequen- 
cies between 1.27 THz and 1.82 THz on both branches 
la and lb. Analyzing the voltage differences AV5 we 
got different values for la and lb (Fig. On branch 
la the voltages of Shapiro steps were detected at regu- 
lar values of hf/ (2e) whereas voltages of lb were shifted 
to values 3% below hf/(2e). This behaviour can be ex- 
plained by assuming that in case lb the junction close 
to the normal electrode becomes resistive while in case 
la a junction inside the stack becomes resistive. This as- 
sumption is also in agreement with the observation that 
the voltage of the second branch at the critical current is 
approximately given by the sum of the voltages of la and 
lb and the voltage differences for the higher branches are 
rather homogenous and are close to the value for branch 
la. 

We want to mention that we could also measure 
Shapiro steps on higher resistive branches. Due to 
slightly varying parameters of different resistive Joseph- 
son junctions analysis of AV s . n is much more complicated 
and was not accurate enough to deduce deviations from 
the second Josephson relation. 

Finally we want to point out that a shift in the Shapiro 
step voltage is only possible, if the resistive junction is 
close to the normal electrode. In recent measurements 
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FIG. 5: IV characteristic of #20. Several current cycles are 
shown in one figure. 
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FIG. 6: Voltages of Shapiro steps measured on IV character- 
istic of sample #20 compared to hf /(2e). The values strongly 
depend on which branch la or lb the resonances occured. 



of Shapiro steps in step-edge junctionals this is different. 
Here the resistive Josephson junction is inside a stack of 
superconducting layers and hence no shift is observed. 
Also no shift of Shapiro steps appears in true 4-point 
measurements and for steps cipssing the zero-current line 
at the point of zero current .t3 



B. Injection experiments in double-mesa structures 

The geometry of samples used for injection experi- 
ments in double-mesa structures is schematically shown 
in Fig. [|. As results obtained for the different samples 
are rather similar, we will discuss only results for sam- 
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FIG. 7: Voltages of Ml and M2 during variation of Ii while 
h = 0. 



FIG. 8: Voltages of Ml and M2 during variation of I\ while 
two junctions of M2 are in the resistive state. 



pie #SR102_3. Here two small mesas Ml of lateral size 
5 x lOfim 2 and M2 of size 4 x lOfim 2 were structured on 
a base mesa B of size 10m x 10/um 2 . The currents 
through the mesas Ml and M2 and the base mesa B and 
the corresponding voltages can be measured sepa- 
rately. The brush-like structure of the IV-curves which is 
similar for the two mesas shows two sets of branches with 
different critical currents. These belong to 10 Josephson 
junctions in Ml and M2 and 6 junctions in the base mesa. 

To explain the operation of the device we want to start 
with the case of no current flowing through M2 (I2 = 0) . 
Measuring V\ and V2 during variation of I\ we obtained 
the curves shown in Fig. |?]. V\(I\) shows the full set of 
IV-curves of Ml and B, while in V<z(J\) only the resis- 
tive junctions of the base mesa B appear. Here V^ii) 
is in fact a 4-point measurement of the IV-characteristic 
of the base mesa. Note that as long as the junctions in 
B are completely superconducting, the voltage V2 is ex- 
actly zero and independent of 1\ . If a small current I2 
is applied to M2, then 1/2(^2) shows the contact resis- 
tance between the gold electrode and the mesa. Again 
the voltage V2 does not depend on I\. 

Now a larger bias current is applied to M2, such that 
some junctions of M2 are in the resistive state. Keeping 
I2 fixed we varied the current I±. During the cycling 
of I\ junctions in Ml are switched on and off into the 
resistive state. From time to time also junctions in the 
other mesa M2 are switched on and off leading to discrete 
jumps in the voltage V2, but otherwise V2 is still constant 
(horizontal line in Fig. |^) . 

In some cases, however, the voltage on M2 jumps to 
an additional branch V^ii), which splits off from the 
constant voltage branch and depends weakly on I\. The 
jump into this branch, which is marked Ay in Fig. ||, 
is always triggered by a switch into one of the higher 
order branches of Ml. Let us note that this happens 
only if I\ and I2 are in opposite direction. In later cycles 
the branch V^/i) can also traced out by increasing the 
current I\ from zero. In general, several sets of voltage 
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FIG. 9: Charge-imbalance generated on the first layer of the 
base mesa by the currents through mesas Ml and M2. 



curves with split branches A V2 (Ii ) occur. The voltages of 
the horizontal branches correspond to different numbers 
of junctions of M2 beeing in the resistive state at the 
fixed current J2. In Fig. ||only the branch corresponding 
to two resistive junctions is shown. 

The maximum voltage difference AV^ is much smaller 
than the voltage between different resistive branches in 
M2 at the bias current I2 ■ It is also much smaller than the 
voltage of resistive junctions in the base mesa. Therefore 
the appearence of AV2(ii) must have a different origin 
and can be explained as follows: If the junctions between 
the lowest layers in both mesas Ml and M2 and the first 
common superconducting layer in B are resistive, then 
a nonequilibrium potential f ^ is generated on the first 
superconducting layer of the base mesa. This is illus- 
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trated in Fig. || In Fig. only the current 1% contributes 
to the charge-imbalance potential, while in Fig. both 
currents contribute. The generated charge-imbalance po- 
tential can be measured directly as additional voltage on 
M2: 



AV 2 (h) = * fl (/i + h) - <M/ 2 ), 



A 2e 2 N{0) 



(17) 



(18) 



Here A is the area of the base mesa. In deriving (|T|) 
we assumed that the charge-imbalance relaxation time is 
large compared to the diffusion time of charge-imbalance 
along the layer. 

Finally we want to explain the asymmetry of 
with respect to the polarity of both currents. The total 
current Ib through B is either 1\ + I2 if the polarity is 
the same, or I\ — 12 if the polarity is opposite. In the first 
case the total current through B is higher. This makes it 
possible that one junction inside B gets resistive before 
the critical current of the small mesas is exceeded. This 
destroys the precondition of a completely superconduct- 
ing mesa B. When the polarity is opposite the current 
through B is always smaller than I\ and junctions in Ml 
will get resistive before any junction in B. 



Determination of the charge-imbalance 
relaxation time 



Both types of experiments can be used to measure the 
charge-imbalance relaxation time. From the shift of the 
Shapiro step we obtain 



SVj2e 2 N(0), 



(19) 



where A is the area of the mesa. From the injection ex- 
periment in the double-mesa structure a similar expres- 
sion is obtained. This equation contains the unknown 
density of states N(0) of the two-dimensional electron gas 
at the Fermi surface. For a rough estimate we may use 
the density of states N(0) — raj (2ith 2 ) of free conduction 
electrons with mass m, then 2e 2 N(0) = 0.67C 2 J _1 m~ 2 . 
Alternatively, by using a — eeo / (2e 2 N (0)d) we can ex- 
press 2e 2 N(0) by the parameter a, the dielectric constant 
e of the barrier and the distance d between superconduct- 
ing layers. In Ref. [l4| the values a = 0.4, e = 20 have 
been estimated from reflectivity experiments (for another 
material). Here it is found that the value of eeo/ (ad) is 
rather close to the value of 2e 2 iV(0) calculated for free 
conduction electrons, which therefore will be used in the 
following. 

With the data of sample #32: A = 64^m 2 , I = 30/^A, 
SV = 50/iV we find r q ~ 70 ps. An estimate for the 



shifts on branch lb of sample #20 with 1.63 THz gives 
values of r q ~ 450 ps. In a similar way we can also 
determine r q from the injection experiments in double- 
mesa structures. Using I\ — 120//A, I2 — 20/iA, AU = 
0.6mV, A = lOO^m 2 we find r q ~ 330 ps, which is similar 
to the values obtained from the shifts of Shapiro steps. 

The charge-imbalance relaxation time r q describes the 
recombination of quasi-particles into Cooper pairs. It 
should be distinguished from the characteristic time for 
the theunalisation of hot quasi-particles, which is much 
shorter .Ej Note that in the present case r q is much 
longer than the period of a Josephson oscillation on a 
Shapiro step. For the charge-imbalance relaxation nor- 
mally inelastic phonon scattering processes are responsi- 
ble. It should be noted that in the case of d-wave pairing 
with an anisotropic gap also elastic impurity scattering 
contributes Jl3 



V. SUMMARY 

In this paper we have discussed new experiments 
showing evidence of non-equilibrium effects in intrinsic 
Josephson contacts in layered superconductors, which 
are due to charge-imbalance produced by a stationary 
bias current. In particular, we have investigated the 
voltage-position of Shapiro steps in the presence of high- 
frequency irradiation. In some cases we observed a down- 
shift of 3% from the canonical value of hf/(2e). This 
shift, which is not a violation of the basic Josephson re- 
lation, can be explained by charge imbalance on the first 
superconducting layer if the Josephson contact is next to 
the normal electrode. 

In another type of experiment we studied the mutual 
influence of currents through two mesas on a common 
base mesa on the measured voltages. The reults can 
be explained by charge-imbalance on the first common 
superconducting layer of the base mesa. This experi- 
ments has .some similarity with the classical experiment 
by Clarkejla where charge imbalance is produced in a 
superconductor by a strong quasi-particle current which 
is then detected as voltage difference between a normal 
contact and a Josephson contact. 

Both experiments allow to measure the charge- 
imbalance relaxation time which is of the order of 100 
ps. 
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